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Abstract. A theoretical study of magnetic field (h) effects on single-particle spectra 
and transport quantities of heavy fermion metals in the paramagnetic phase is carried 
out. We have employed a non-perturbative local moment approach (LMA) to the 
asymmetric periodic Anderson model within the dynamical mean field framework. 
The lattice coherence scale w^, which is proportional within the LMA to the spin- flip 
energy scale, and has been shown in earlier studies to be the energy scale at which 
crossover to single impurity physics occurs, increases monotonically with increasing 
magnetic field. The many body Kondo resonance in the density of states at the 
Fermi level splits into two with the splitting being proportional to the field itself. 
For h> 0, we demonstrate adiabatic continuity from the strongly interacting case to 
a corresponding non- interacting limit, thus establishing Fermi liquid behaviour for 
heavy fermion metals in the presence of magnetic field. In the Kondo lattice regime, 
the theoretically computed magnetoresistance is found to be negative in the entire 
temperature range. We argue that such a result could be understood at T > cj^ 
by field- induced suppression of spin- flip scattering and at T < through lattice 
coherence. The coherence peak in the heavy fermion resistivity diminishes and moves 
to higher temperatures with increasing field. Direct comparison of the theoretical 
results to the field dependent resistivity measurements in CeBg yields good agreement. 

PACS numbers: 71.27.-Fa,71.28.-Fd,71.30.-Fh,75.20.Hr 
1. Introduction 

The investigation of lanthanide/actinide based heavy fermion(HF) systems has been 
a central theme in condensed matter physics both theoretically and experimentally [1]. 
Their behaviour is quite distinct from conventional clean metals, the basic physics being 
driven by strong spin-flip scattering from essentially localized f-levels, generating the 
large effective mass. The periodic Anderson model(PAM) forms the general paradigm 
within which these materials are studied. The minimal model consists of a regular array 
of sites, each associated with a localized, non-degenerate f-electron core orbital, coupled 
to a delocalized conduction electron orbital via a local hybridization. Neighboring 
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conduction electron orbitals are connected via a hopping matrix element and electron 
interactions enter the model via an on-site coulomb repulsion. 

The dynamical mean field theory (DMFT) [2] has proved to be a very powerful 
framework for studies of various lattice models such as the Hubbard model or the PAM. 
Within DMFT, which is exact in the limit of infinite dimensions, the self-energy becomes 
spatially local or momentum independent. As a consequence, lattice models map onto 
an effective single-impurity Anderson model (SIAM) with a self consistently determined 
host [2]. 

The PAM in the absence of a magnetic field has been studied extensively within 
the framework of DMFT. Some of the numerical or semi- analytical methods (impurity 
solvers) that exist for solving the the effective SIAM that arises within DMFT are 
quantum Monte Carlo [3], numerical renormalization group [4], iterated perturbation 
theory [5], noncrossing approximation [6], the local moment approach [7,8], the large- 
N/slave-boson [9], exact diagonalization [10] self-consistent perturbation theory [11], the 
Gutzwiller variational method [12] and the average T- matrix apprximation [13]. Every 
method has its own advantage, however most of them suffer from one or the other 
limitations. For example, the QMC suffers from minus sign problem at low temperature 
or large interactions, while IPT is able to capture only an algebraic decay of the Kondo 
scale with increasing interactions. ED and NRG are in principle exact, as is QMC, 
but since the spectral functions arc obtained as a set of discrete poles, a broadening is 
required, which is non-uniqucly specified. 

In this context, the (LMA) [8, 14-18], a diagrammatic theory based non- 
perturbative many body method, has emerged as an approach for the single impurity 
model and even for the lattice models within DMFT, that overcomes some of the 
limitations mentioned above. In particular, one obtains the exact dependence on 
interactions and hybridization of the Kondo scale of the SIAM in strong couphng as 
the Bethe ansatz solution for the Kondo model [14]. Very good agreement has been 
seen between LMA and NRG results for the spectral functions of the SIAM [17, 18]. 
Within DMFT, the LMA has proved quite successful in describing the spectral and 
transport properties of several metallic and insulating heavy fermion systems such as 
CeBg, SmBg, YbAls, CeAls, YbBi2, CeOs4Sbi2 in the paramagnetic phase [19-21]. The 
LMA is computationally inexpensive, yields quantities in real frequency directly, and is 
semi-analytical. One limitation of LMA is that it is based on a symmetry restoration 
ansatz, that is not easily generalizable to other problems, such as the multi-orbital or 
multi-channel cases. Recently, a generalization of LMA has been reported for the multi- 
orbital Anderson and Hubbard model [22]. Further, it was shown that the LMA is a 
conserving approximation. The other limitation is that since LMA is an approximate 
theory, one has to benchmark its results against more exact theories such as NRG and 
QMC to ascertain its reliability. Nevertheless, given the several advantages above, and 
benchmarks, the LMA within DMFT is an appropriate choice to study the effects of 
magnetic fields in heavy fermions. 

Several theoretical studies of the effects of magnetic field on heavy fermion systems 
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using either the Kondo lattice model (KLM) or the PAM have been reported. The 
magnetic field induced insulator-metal transition in Kondo insulators has been studied 
by various groups [23-25]. As relevant to heavy fermion metals, the metamagnetic 
transition has been studied using large-N mean-field [26] and subsequently using 
DMFT-I-QMC [27]. In this work, our objective is to understand the magnetotransport 
in heavy fermion metals. Previous work in this direction has been carried out mainly 
either using the single-site Anderson models [28] , thus missing out the lattice coherence 
effects completely or using large-N mean field treatment of the Anderson lattice 
Hamiltonian [29]. In order to capture the lattice coherence effects along with the 
single-impurity incoherent regime quantitatively within a single framework, wc employ 
the finite-field LMA [17, 18, 25] within DMFT for the periodic Anderson model away 
from half-filling and determine the effect of magnetic field on spectra and transport 
of heavy fermion metals. Our focus has been on the strong coupling regime, and the 
quantities that we have studied are spectral functions and magnetoresistance. The paper 
is organized as follows; the model and formalism are presented in section II, followed, in 
section III by results for the field evolution of single-particle dynamics and d.c transport. 
In section IV, experimental magnetoresistance of CbBq is compared with theory; and 
finally we conclude with a brief summary. 

2. Model and formalism 

The Hamiltonian for the PAM in standard notation is given by: 

^ = E 4s-. + E(^/ + ^flJi-MLu + V Y.(fl^i^ + ^■^) + E ^c4c,.(2.i) 

{i,j),u icr ia ia 

The first term describes the kinetic energy of the nonintcracting conduction (c) band 
due to nearest neighbour hopping t. The second term refers to the /-levels with site 
energies e/ and on-site repulsion U, while the third term describes the c// hybridization 
via the local matrix element V. The final term represents the c-electron orbital energy. 
In the limit of large dimensions, the hopping needs to be scaled ast oc t*/ \/Z, where Z is 
the lattice coordination number. We consider the hypercubic lattice, for-which the non- 
interacting density of states is an unbounded Gaussian {po{e) — exp{—e^/t'^^/-\/~{'Kt^)). 
Particle-hole asymmetry in the PAM could be introduced in two ways [8]: (i) through 
an asymmetric conduction band, i.e. ec 7^ 0, or (ii) through an asymmetric f-level, 
e/ 7^ — In general, we may quantify the f-level asymmetry by defining r] = 1 + 
such that ?7 = is equivalent to particle-hole symmetric f-levels. Our primary interest 
is in the strong coupling Kondo lattice regime (n/ — > 1) but with arbitrary conduction 
band filling (ric). 

Within DMFT, the PAM may be mapped onto an effective self-consistent impurity 
problem within DMFT. We choose the local moment approach to solve the effective 
impurity problem arising within DMFT. For details of LMA developed for use within 
DMFT in the absence of a field, the reader is referred to some of our previous 
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work [8,19-21]. Magnetic field effects in the symmetric PAM as appropriate to Kondo 
insulators was recently studied by us [25]. As discussed in that work, the presence of a 
global magnetic field results in the Zeeman splitting of the bare electronic energy levels 
as e^o- = f 7 ~ crhy, for 7 = c and / electrons. Here = ^g^fisH and /is is the Bohr 
magneton; the constants gj and gc are the electronic g- factors for the / and c electrons 
respectively. Although gf gc in general, for simplicity we set gf = gc- The degeneracy 
of the symmetry broken solutions at the mean-field level (denoted by A and B) is lifted 
by a magnetic field. We consider here /i > for which (A-type) is the sole 

solution. The Feenberg self energy [8, 19, 25] becomes a functional solely of the A-type 
Green's function, i.e. Sa{uj) = Sa[G'^^^]. Since the B-type solution does not exist, the 
label 'A' will be implicit in the following. 

In presence of a uniform magnetic field, the LMA Green functions are given by 

,+ , ah-ec- 5^[G^] - 



Gi{uj,T,h) 



uj^ + ah — Sf — Eo.(a;, T, h) — 



a;+ + ah — €f — ^^.(a;, T, h) 



(2.2) 



(2.3) 



with G'^{ijj) — I Ylia ^Zi^)- The self energy can be written as a sum of the static and 
dynamic parts as. 



Ecr(a;,T, h) = ^(n - a\fi\{T, h)) + E^{uj,T, h) 



(2.4) 



where |/u|(T, h) is the UHF local moment in presence of magnetic field and temperature. 
The dynamical self-energy within the LMA is given by [16,25] 

V„{uj2, h) 



/OO /"OO 
/ duJr(Lj2x""'{oJi,T,h) 
-00 J -00 

/i(a;i; 0^2) = ^^(-0^1) [1 - f{u:2, T)] + e{ujr)f{uj2, T) 



h{uJx]UJ2) 



(2.5) 



where x~^"{'^^T^ h) = n'Hmll^'^'^ (uj , h), f(uj) = [exp'^'^ +1]~^ is the Fermi function and 
9{uj) is the Hcavisidc step-function. Here, Il~'^'^{uj,T, h) denotes the transverse spin 
polarization propagator which can be expressed as 11'*'" = °n''"~/(l — U^~^~) where 
°Q''"~, the bare p-h bubble, is constructed using the field dependent mean-field spectral 
densities [14,20]. The host spectral function is given by 'D„{u,h) = — ^Im^o-(c<;, /i); 
where the host/medium Green's function Q„ is given by, 

r2 H-l 



Qa{oJ,T, h) = 



Cf + ax + ah 



(2.6) 



u;+ -Sc + ah- S^iu, T, h) _ 

where the parameters x — U\ix\/2 and e/ are determined at /i,!" = by satisfying the 
symmetry restoration condition 

Sf (cu = 0; e/, x) - Sf (cj = 0; e/, x) = t/|/x(e/, x)\. (2.7) 

and the Luttinger's integral theorem [8,30] 

du (9E(u;) 



h{ef,x)^lmJ_^-^Gf(u) 







(2.8) 
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which in turn ensures a Fermi hquid ground state. Note that ej is distinct from Sf. 
The latter is a 'bare' model parameter while the former is a derived 'shifted chemical 
potential', which is determined through the imposition of the Luttinger's integral, 
or equivalently the Friedel sum rule [8]. In practice however, it is more convenient 
numerically to fix x and e/ at the outset and treat U and e/ as unknown parameters 
to be determined by the above two conditions [8]. As input for the calculation at 
a given field/temperature, we use the self-energies and Green's functions of a lower 
field / temperature. 

In the limit of infinite dimensions, vertex corrections in the skeleton expansion for 

the current-current correlation function are absent, hence a knowledge of single-particle 
dynamics is sufficient within DMFT to determine q=0 transport properties. For h > 0, 
the (i-dimensional isotropic conductivity is computed by adding the contributions from 
each of the spin channels as 

aiu, T,h)^^J2 du;^[f{u;^) - f{u;, + cu)] (L'^(e; u;^)Dl{e■, cu, + cu)), (2.9) 

where o"o = ^ (a is the lattice parameters) and D^(e;uj) = — ^ImG'^(e; a;). The lattice 
c- Green's function is given by G^(e;a;) = [ja-{^^,T, h) — e]~^ where 

-f^^oo^ + ah-ec- (ou'^ + ah - Sf - t^{ou, T, h)^ ^ , 
and the e - average is defined by, 

/oo 
depo{e)A{e;u;) (2.10) 
-oo 

We now proceed to discuss the results obtained by implementing the finite field and 
finite temperature LMA. 

3. Results and discussions 

In this section, we will discuss the single particle dynamics and transport in presence 
of magnetic field using LMA. Our primary focus is on the strong coupling Kondo 
lattice regime (where Uf 1, but Uc is arbitrary). Before discussing the finite field 
calculations, we will review a few well established concepts. At T = and h = [8], 
the strong couphng Kondo lattice regime is characterized by an exponentially small (in 
strong coupling) low- energy scale lol — ZV'^ jt^. The single-particle properties of the 
asymmetric PAM exhibit universal scaling in terms oiuj/ujL and for a fixed Cc/t^ and 77. 
For finite-T and h = Q, the spectra D'^{uj) and V^D^{uj) exhibit scaling in terms oiuj /ujl 
and T/Tl{Tl = ujl). In summary, the spectra and transport properties of the PAM are 
universal functions of uj/uj^ and T/ujl, for a given conduction band filling and f-level 
asymmetry; thus being independent of the bare f//t* and V/t*. In this context, it is 
important to mention that the universal form of the scaling functions does depend, 
albeit weakly, on the specific lattice, which manifests in the bare conduction band 
density of states. The apphcation of magnetic field does not destroy this universality as 
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shown in a recent work on the symmetric PAM (T = 0) where such scahng in terms of 
Uh = ujt^:/{Z{h)V'^) for a fixed effective field (detailed description is in [25]) was shown 
to hold good in the presence of field as well. 

To see if such scaling occurs for the asymmetric case, we carry out a low frequency 
Fermi liquid analysis of the Green's functions by expanding the self-energy about the 
Fermi level to first order in u as, 

S«(.;, h) = S^(0, h) + {l- ^)u; (3.1) 

Substituting equation (13. ip in equations (12.21) and (12. 3p . we find that the spin dependent 
spectral functions are just renormalized versions of their non-interacting counterparts 
and are given by(neglecting the 'bare' terms h and u in the strong coupling limit), 

D^u; h) --4° po I -6. - ^ : 1 (3.2) 



V\ijJha - e} + ^Kff) 

where ujho- = ^t^,/ {Zcj{h)V'^) and the renormalized f-level (e^) and hlj^ can be expressed 



as, 



^f = Y72i^f + ^ -^ + S^(0,0)) (3.4) 



and 

h"/^ = y-2ih- a(Sf (0, h) - S^(0, 0))) (3.5) 

Using the SR condition, equation (12. 7p . it is easy to see that the effective /-level, e*j is 
spin independent. Note that the quantities thus defined, namely, Uha, and h*^^^ are 
dimensionless. 

From equations (13. 2p and (13. 3p . we can infer the following: (i) The spin dependent 
spectra D'^{ijj,h) and V'^Dl{uj,h) should exhibit scaling in terms of uJha for a fixed 
h^Jf i.e, for a fixed h'^-^^, if we plot the spin dependent spectra vs uJha for different 
values of U, they should collapse onto a single curve at the low energy regions for a 
particular h^^^. This is shown in figure [1] where in the insets, the conduction electron 
spectra for up spin (top panel) and the f-electron spectra for up spin (bottom panel) 
are plotted as a function of the bare frequency u/t^ with ec = 0.3t*, i] ^ 0, V'^ = 0.2tl 
and T=0 for three parameter sets: U ~ 5.1^=,, and h = 0.3 (solid), U ~ 6.1^=,, and 
h = 0.25 (dotted), U ~ 6.6t„ and h = 0.24 (dashed) such that /if ^=0.39 is fixed (where 
h = ht^,/{Z{0)V'^)). In the main panel the same spectra are plotted as a function of 
LJhi- We see that when we plot the spectra vs bare frequency, they appear very different, 
but when plotted vs. ujhi (main panel: top and bottom), they collapse onto a single 
universal form. We have checked that similar universality holds for the down spin also 
for a fixed h^/'^ . 

(ii) The spin dependent conduction and f-electron spectra D^{uj) and V'^D^u) should 
adiabatically connect to the non-interacting limit at low energy scales i.e, if we take 
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Figure 1. Insets: LMA spin dependent conduction and f- electron spectra vs bare 
frequency uj/t^ for Cc — 0.3t*, rj ^ 0, = 0.2tl and T = for three parameter sets: 
U ~ 5.1^, and h = 0.3 (solid), U - 6.1i, and h = 0.25 (dotted), U ~ 6.6^,, and 
h = 0.24 (dashed) for fixed hf^' = 0.39 ( where h = ht^/{Z{0)V'^)). Main: The same 
spectra when plotted vs ujhi — ojt^l {Z^(h)V'^ collapse into one single universal form. 

the c-electron and /-electron fields as zero and h^J^ respectively, substitute for e/ the 
renormalized f-level (e^) and compute the spectra in the non-interacting limit; then 
the interacting spectra and the non-interacting spectra should be same at low energy 
scales. This is shown in figure |21 where the spin dependent scaling spectra D'^{ijj) for the 
interacting case {U = 6.6t*) is superposed onto the non-interacting spectra. We see from 
figure[2]that both the curves are almost identical near the Fermi level. This demonstrates 
adiabatic continuity of the strong coupling regime to the non interacting limit which 
represents Fermi liquid behaviour in presence of magnetic field for the asymmetric case 
(in parallel to the symmetric PAM [25]). 

To see the effect of magnetic field on the spin dependent effective fields h"^^^ and the 
spin dependent quasiparticle weights Z^j, we have plotted in figure [3], the spin dependent 
effective fields [h^^ and h'j^^) vs h (left panel) and the spin dependent quasiparticle 
weights {Z^{h) and Z^{h)) vs h (right panel). We see that Z^{h) = Z^{h) and h^"^ = 
up to about h ~ 10.0. So, we can assume that the spin summed spectra D^{ijj) and 
V'^D^iuj) also should exhibit scaling in terms of Uha for low fields (up to /i ~ 10) 
for a fixed h'^J^ . To see this, we plot up spin scaling spectra (main panel) as well as 
corresponding spin summed spectra in figure H] (insets) for U = 5.1t^ (solid line) and 
U = 6.6t* (dotted line) for a fixed hf^ = 4.2 {h ~ 10.0) and T = 0. Indeed, we see that 
the spin summed spectra also exhibit universal scaling and adiabatic continuity at low 
energy scales. We wish to reiterate the dicussion in the beginning of this section that 
that the universality and scaling demonstrated above holds for a fixed conduction band 
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Figure 2. LMA spin dependent scaling spectra (solid lines) are superposed onto 
the corresponding non-interacting spectra (dotted lines) for a fixed h^'^ = 0-45 and 
e*f = 1.09. 




h h 

Figure 3. Left panel: h'^^^ vsh — ht^/ {ZV'^){ Z is the quasiparticle weight for h = 0) 
for Cc = 0.3i*, rj 0^ T = and U = 6.6t». Right panel: Spin dependent quasiparticle 
weight Za-{h) vs h for the same parameters. 

Uc, a fixed /-level asymmetry and for a specific bare conduction band density of 
states. In other words, the functional forms of the spectra shown above would change 
for a different choice for or t] or po{e). 

Now we turn to the LMA results for the field evolution of the asymmetric PAM 
for a fixed temperature. In figure [5], we show T = spin summed conduction and /- 
electron scaling spectra for various fields: h = (solid), 0.5 (dotted), 1.0 (dashed), 1.5 
(dot-dashed), 2.0 (bold sohd), 3.0 (bold dotted), 5.0 (bold dashed) and 10.0 (bold dot- 
dashed) and fixed interaction strength U = 6.6t*. In our earlier work on the symmetric 
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Figure 4. Main panel: LMA up-spin conduction electron scaling spectra (top) and up 
spin /-electron scaling spectra (bottom) vs ojhi for Cc = 0.3t*, 77 — > and V"^ = 0.2tl 
for two parameter sets: U ~ 5.1^* and h = 10.0 (solid), U ~ 6.6t, and h = 9.0 (dotted) 
for fixed h'^^^ = 4.2. Insets: The corresponding spin summed conduction and (top) /- 
electron (bottom) scaling spectra vs ujh] for the same parameter sets. Spin summed as 
well as the spin dependent spectra exhibit universality/scaling in the strong coupling 
regime. 

PAM [25], we have shown that there is a gap at the Fermi level, as there must be for a 
Kondo insulator. But in the asymmetric limit ec = QM^.h = (see figure [5]), the gap 
moves away from the Fermi level and becomes a pseudo-gap. The width of the Kondo 
resonance at the Fermi level is proportional to the low-energy scale ujl„ = Z„{0)V'^. 
The insets of figure [5] show the spectra on a large frequency scale. It is seen that 
for very large frequencies the tails of the spectra for all fields are identical which is 
physically natural since one can expect that the effect of the field should dominate 
only for < ht^/{Z^{h)V^). Now with increasing magnetic field, the lattice Kondo 
resonance splits into two peaks, with the distance between the peaks also increasing. 
Qualitatively we can understand the results as follows: In the non-interacting limit for 
the symmetric case, we have shown that the c and / -levels shift rigidly due to the 
Zeeman effect. The same concept is valid for the asymmetric case also for the non- 
interacting limit i.e. the spectral function for the up and down spin bands shift rigidly 
away from the Fermi level. In the presence of interactions, in parallel to the symmetric 
case, the shift of the spin bands should not be rigid due to the competition between 
Zeeman splitting and Kondo screening. And indeed, we find that although the distance 
between the two peaks varies linearly with field in strong coupling (see figure [6]) , the 
slope is not equal to 2 as would have been if the shift had been rigid (or purely due 
to Zeeman effect). In fact, the slope should be equal to 4 as the following argument 
shows. From equations (13. 2p . (I3.3p . (13.41) and (13. 5p . it is straightforward to see that 
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Figure 5. LMA conduction and /- electron spectra vs Cjh-^ — Lot^,/{Z^{h)V'^) for 
U = 6.6i,, ec = 0.3i, and 77 ^ 0.0 for = (solid), h = 0.5 (dotted), 1.0 (dashed), 
1.5 (dot-dashed), 2.0 (bold solid), 3.0 (bold dotted), 5.0 (bold dashed) and 10.0 (bold 
dot-dashed) for T = 0. 



the splitting of the Kondo resonance, Skr is given by 

Skr = Za{h) [h - a (S^(0, h) - S^(0, 0))] (3.6) 

cr 

The strong-couphng asymptotic behaviour of the above quantity may be inferred using 
the expressions derived in [17] for the flat-band SIAM. 



S^^(0,/.)-S,«(0,0) = -^ln 

vr 



(3.7) 



Zcr{h) 71 h 

Z^M^^^ 2AoZ.(0) ^^-^^ 

where Aq = ttV'^Pq{0) is the flat-band hybridization. Although these were derived for 
the SIAM, we are justified in using the same expressions for the PAM for the following 
reason: Since the transverse spin-polarization propagator is constructed using the UHF 
propagators which do not contain the low energy scale, and hence are flat for ~ C^(l), 
the strong coupling (SC) asymptotics derived for the SIAM would be similar to that 
for the PAM. The only difference would be in the value of Aq ~ 0{V'^/t^,) since the 
UHF propagators for the PAM have a different structure than that of the SIAM. Using 
equations (13.71) and (I3.8P in (13.61) . and expanding the logarithm appearing in equation 
( 13. Sp to linear order in h, we get 

Skr ^ 4/1 (3.9) 
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As argued in [17], the above result is synonymous to the Wilson's ratio being equal 
to 2 in SC. The inset of figure [6] shows the ratio SxR/h as a function of interaction 
strength. It is seen that even for ?7 ~ 2, the ratio is ~ 3, which implies a non-rigid 
shift of the spin-bands. Further, with increasing interaction strength, the ratio increases 
monotonically. Although it must asymptotically approach 4, to be consistent with the 
result (equation (13. 9p ) obtained above, we are unable to access the large U region (for 
ec = 0.5) due to the prohibitive computational expense in handling exponentially small 
low energy scales. 




Figure 6. Splitting of the Kondo resonance vs field for U = 5.5t, and Cc — 0.5, V"^ = 
0.2tJ. The slope of the line in the main panel, i.e. Skb/^ a function of interaction 
strength U is shown in the inset. 

In a recent DMFT study using QMC as the impurity solver, the Kondo lattice 
model was studied [27]. Real frequency spectra were obtained with a stochastic analytic 
continuation method. The parameters chosen were J = 1.6t* and Uc = 0.85. The 
main result was that the bands were found to shift rigidly, in contrast to what we 
find above. The difference with our findings could be due to various reasons, of which 
the most important seems to be that their low energy scale is 0.09t^,, while our scales 
are ~ 10~^t* (see figure 3). This is significant because only when the Kondo scale 
is exponentially small is the coupling between the impurity and the conduction spin 
renormalized strongly, and thus the Kondo screening would be strong. If this Kondo 
screening is weak, as could probably be in the study mentioned above, then the Zeeman 
effect would win easily, and the bands would shift rigidly. 

Now, we turn our attention to field-dependent transport properties. At finite 
temperature, as the spin summed conduction electron spectra exhibits universal scaling 
in terms of T/uia {j^Lcr = Z„{h)V'^ /t^,) in the strong coupling regime for low fields, 
we expect that the resistivity will also exhibit scaling in terms of T/uia- The classic 
HF metallic resistivity increases with temperatures (initially as T^), goes through a 
maximum Tmax — (the peak position of p(T)) and then decreases with temperature 
characteristic of single impurity behaviour. We observe the same form for the resistivity 
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Figure 7. Top: Left panel: Resistivity vs ojhi for the same parameters, but for the 
various fields {h = 0.0 (solid), 0.06 (dotted), 0.1 (dashed), 0.5 (long dashed) and 2.0 
(dot-dashed). Right panel: Resistivity ws for the same sets of parameters. Bottom: 
Magnetoresistivity vs bare fields /i/t* for U = 6.6t,, Cc = 0.3t* and 77 — > for T = 0.3 
(sohd), 0.6 (dotted), 2.0 (dashed), 6.0 (dot-dashed). 

shown for both a =t, i and for a range of fields: h = 0.0 (solid), 0.06 (dotted), 0.5 
(dashed), 0.5 (long-dashed) and 2.0 (dot-dashed) with U/t^ equal to 6.6. The two 
most significant features of the above result are - (i) Negative magnetoresistance is 
observed over the entire range of temperatures as shown in figure [7] (bottom), (ii) 
The maximum in resistivity at the coherence peak reduces in magnitude and moves 
to higher temperatures with increasing field. In previous work, magnetoresistance 
was found to be positive for the SIAM [31] through Fermi liquid theory arguments 
and for the Anderson lattice [13] through the average T- matrix approximation (ATA). 
The fundamental difference between the SIAM and the PAM is the presence of lattice 
coherence in the latter which manifests itself at low temperatures. At T > ul, it is 
clear that magnetic field suppresses spin-flip scattering (SFS) since the main effect of 
a field would be to polarize the system. However, at low temperatures, the physics 
is different. Kondo singlet formation in the SIAM, which at zero field, quenches the 
SFS and hence leads to saturation of resistivity, is inhibited at finite fields, implying an 
increase in incoherent scattering and thus a positive magnetoresistance. In the lattice 
case, however, although magnetic field does inhibit the singlet-forming screening process 
of local moments, the periodicity of the local moments (even though not fully screened), 
introduces lattice coherence which strongly suppresses incoherent scattering at low 
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temperatures, and hence again leads to negative magnetoresistance. In other words, the 
incoherent scattering introduced by magnetic field in the SIAM through the inhibition of 
the screening process, is countered in the PAM because of the presence of lattice/Bloch 
coherence. It is possible that such an effect is enhanced in the present calculation due to 
the use of the DMFT framework, since non-local dynamical fiuctuations are completely 
neglected in this framework. The ATA for the Anderson lattice is in a subtle way 
somewhat similar to the the single-site DMFT without the self- consistency. Which would 
imply that the effects of lattice coherence are probably suppressed in the ATA, and hence 
positive magnetoresistance is observed. Thus, we conclude that, at higher temperatures 
(T > ujl), suppression of the spin-flip scattering, and at lower temperatures, the presence 
of lattice coherence in the PAM lead to negative magnetoresistance in the presence of a 
field for all T. 

The shift of coherence peak to higher T with increasing field is a refiection of the 
increase in the low energy scale with field (see figure 3). The latter result has been 
obtained for the SIAM as well [17]. Since increasing interactions lead to an exponential 
decrease in uj^, it appears that the polarizing effect of magnetic field is to counter the 
effect of interactions, and eventually at very large fields, wipe out local moment physics 
(scattering/screening) completely. At temperatures much higher than the field(T ^ h), 
the effect of magnetic field is negligible implying that the magnetoresistance is almost 
zero. Next, we compare our theoretical results with experiment. 

4. Comparison to CcBq 

In this section we want to compare our theoretical results with experiment on CcBq. 
The rare-earth hexaboride CcBq has been investigated for many years [32-37]. The 
cubic-lattice system, at low temperatures, exhibits various magnetic phase transitions 
between 1.6-3.3K, which manifest clearly as kinks in the resistivity, above which the 
system is in a paramagnetic phase. It is the T >3.3K phase that we concentrate on, 
since our approach does not describe the symmetry broken states. Although there have 
been extensive studies of this material, the most detailed magnetoresistance study was 
carried out by Takase et al. [32], who measured the resistivity of CcBq single crystal in 
the temperature range from 3 to 300K and up to magnetic field 85 KOe. We use their 
data for comparison to our theory. 

From an earlier study [21], it is known that the HF system CcBq, belongs to 
the moderately strong coupling regime. Hence we have taken the theoretical results 
to compare with experiments on CcBq for the following parameters: U = 2.4t^,, 
Cc = est*, ?7 and V'^ = 0.2tl. The experimental resistivity for zero field shown 
in the left panel of figure [H](data from Takase et al. [32]) is characteristic of classic 
HF metals. The resistivity rises sharply from a low value, going through a coherence 
peak at T ~4K, and subsequently decreasing through a small log-linear regime. At 
higher temperatures [21, 33] (not shown), resistivity exhibits single-impurity incoherent 
behaviour (T ~ 100 — 300K), goes through a weak minimum (T ~ 375K) and finally 
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starts increasing again like a normal metal. With increasing magnetic field, negative 
magnetoresistance is seen at low temperatures (T < 50K) for all fields. At higher 
temperatures, the resistivity remains unaffected by magnetic field. This behaviour 
is natural and expected as discussed in section 3. The application of magnetic field 
results in suppression of spin-spin scattering at T > ul, thus causing a reduction in 
resistivity. At low temperatures lattice coherence takes over and the magnetoresistance 
stays negative. In the strong coupling regime, the effect of a magnetic field (H) should be 
expected to extend to temperatures, which are of a similar magnitude i.e T ~ gfiBH/ks, 
which in the experiment corresponds to ~ 5.6K for H = 85kOe (see next paragraph). 
However, the highest field affects the resistivity upto a temperature of ~ 50K. Such 
behaviour is characteristic of intermediate coupling regime [20], which is consistent with 
the model parameters for CeBg. 

To compare the experimental results with our theory, two fundamental requirements 
need to be met. First requirement is that, we should extract the contribution to the 
measured resistivity from phonons(pp/i(T)) and the residual resistivity(p(0)). This is 
given by, p^^gi'^) — '^(p(^) ~ P(0)) ~ PphiT) where, a is a constant which comes from 
the error of the sample geometry to the measured resistivity. A detailed discussion of 
this point is given in ref [21]. In this work, we simply assume a = 1. The second 
is, the value of low-energy scale is needed for comparing our theoretical results with 
experiment. For this, we superposed the theoretical resistivity pmagiT) onto PmagC^) ^'^^ 
h = to calculate the value of the low-energy scale. The value of the low energy scale 
turns out to be ul = 2.2K. Assuming the (^-factor to be roughly unity, the magnetic 
field of IkOe can be translated into multiples of the low energy scale. For H = IkOe, 
\gpBH ~ 0.033/^ ~ 0.015a;L. Thus the magnetic fields employed in the experiments [32] 
turn out to be (in multiples of ojl), h = 0, 0.3, 0.6, 0.9 and 1.25 corresponding to fields of 
H = 0, 20kOe, 40kOe, 60kOe and 85kOe respectively. So, given the model parameters 
and the low energy scale, along with the temperature range and the field range, we 
can compute the resistivity as a function of temperature at the same fields as in the 
experiment. 

With a simple multiplicative scaling of the theoretical resistivities on the x-axis by 
the low energy scale and the y-axis by a single multiplicative factor (same for all fields), 
we show the theoretical computed dc resistivities for the same magnetic fields as the 
experiment in the right panel of figure [H The experimental resistivity is found to rise 
more steeply than the theoretical one. The magnitude of magnetoresistivity is found 
to be higher in theory than the experiment. Nevertheless the functional form agrees 
quite well. At T > 50K, the magnetic field is seen to have very little effect on the 
resistivity, both in theory and experiment. Below, 50K, the magnetoresistance in the 
theory is negative, which is also in agreement with experiment. The coherence peak 
is seen to move to higher T with an increase in field, which behaviour is seen in the 
experiment as well. The small log-linear regime that appears at temperatures higher 
than the coherence peak temperature is also reproduced in the theory. Thus, when 
compared with the experimental results in figure [H good qualitative agreement is found 




Figure 8. Resistivity vs temperature of CeB^ (a) experiment (b) theory for the 
various fields: H = Q kOe (sohd), 20 kOe (dotted), 40 kOe (dashed), 60 kOe (long 
dashed) and 85 kOe (dot dashed). 



between theory and experiment. 



5. Conclusions. 



In summary, we have employed a non-pcrturbative local moment approach to the 
asymmetric periodic Anderson model within DMFT in presence of magnetic field. Field- 
dependent dynamics and transport properties of the model have been computed. In the 
strong coupling Kondo lattice regime of the model, the local c- and /-electron spectral 
functions are found to exhibit universal scaling, being functions solely of uj/ul^, T/ujl^ 
{ujLa being the low-energy scale) for a given effective field /ig//- Although the externally 
applied field is globally uniform, the effective local field experienced by the c- and 
/-electrons differs because of correlation effects. Fermi hquid behaviour has been 
established even in presence of magnetic fields through adiabatic continuity to the 
non-interacting limit. In presence of magnetic field, the quasiparticle peak at /i = 
and T = splits into two. The shift of these peaks away from the Fermi level is 
not rigid due to the competition between local moment screening and Zeeman spin- 
polarization. Although these shifts vary linearly with the field in strong coupling, the 
slope is enhanced as compared to the the non-interacting limit. Finally, a comparison 
of theoretical magnetoresistance results with those of CeBg, measured experimentally, 
yields good agreement. 
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